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Abstract 


The present study focuses on the basal clade of the genus Poecilimon — subgenus 
Hamatopoecilimon, with a complex distribution in the eastern Aegean. We com- 
pare different mitochondrial gene trees, test the phylogenetic utility of the mito- 
chondrial A+T-rich region and discuss mitochondrial phylogeny. We infer that the 
ancestral lineages of Hamatopoecilimon split rather simultaneously, affected by 
geological events during late Miocene. 


Zusammenfassung 


Die vorliegende Arbeit konzentriert sich auf die basale Gruppe der Gattung Poe- 
cilimon mit der Untergattung Hamatopoecilimon, die eine komplexe Verbreitung in 
der Ostagais aufweist. Wir vergleichen verschiedene Stammbaume, basierend auf 
Sequenzen mitochondrialer Gene, testen die mitochondriale A-T-reiche Region 
auf phylogenetische Aussagekraft und diskutieren die mitochondriale Phylogenie 
der Gruppe. Wir nehmen an, dass die ursprünglichen Linien von Hamatopoecili- 
mon nahezu gleichzeitig entstanden sind, beeinflusst durch die geologischen Vor- 
gange im spaten Miozan. 


Introduction 


Poecilimon, Fischer 1853 is the most diverse bush-cricket genus in the Palearctic 
including more than 140 species (CIGLIANO et al. 2020). Despite the wide range of 
several species, diversity is concentrated in the Balkans, Asia Minor and the Cau- 
casus where levels of endemism are noticeably high. This pattern is explained by 
rapid diversification, guided by geological events in the region during the Neogene 
and by the following Pleistocene climatic cycles (HEWITT 1996, LA GRECA 1999, 
ÇIPLAK 2004, 2010, KAYA 2015). 


Variety of broad studies accumulated knowledge about the distribution, morphol- 
ogy, bioacoustics, karyology and systematics of the genus (RAMME 1933, BEY- 
BIENKO 1954, HELLER 1984, 1988, WARCHALOWSKA-SLIVA et al. 2000, ULLRICH et 
al. 2010, UNAL 2010, GRZYWACS et al. 2014). Recently authors focus their research 
on separate species groups which share common morpho-acoustic traits and evo- 
lutionary history (HELLER et al. 2005, 2006, 2011, HELLER & SEVGILI 2005, CHO- 
BANOV & HELLER 2010, KAYA 2018, KAYA et al. 2012, 2015, 2018, SEVGILI et al. 
2018). 


A phylogenetic study (ULLRICH et al. 2010) based on mitochondrial and nuclear 
markers suggested update to the relationships within Poecilimon with Poecilimon 
hamatus species group representing its basal branch. 
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A detailed research on morphology and acoustic behavior in the group followed 
(HELLER et al. 2011), raising the group rank to a subgenus level — Hamatopoe- 
cilimon. Six species are recognized within the subgenus — P. deplanatus Brunner 
von Wattenwyl, 1891, P. hamatus Brunner von Wattenwyl, 1878, P. ikariensis Wil- 
lemse, 1982, P. klausgerhardi Fontana, 2004, P. paros Heller & Reinhold, 1992, 
and P. unispinosus Brunner von Wattenwyl, 1878. All are found over Aegean is- 
lands and the mainland of western Anatolia with three species endemic to single 
islands (See HELLER et al. 2011 for detailed distribution map). Distribution of the 
subgenus has been thus related to the paleogeographic evolution of the region — 
the Neogene disintegration of the Aegeid plate, which triggers speciation by vicar- 
iance, and the subsequent reconnection of landmasses, which could reestablish 
contact between populations (CIPLAK et al. 2010, POULAKAKIS et al. 2014). 


Our study aims to reconstruct a robust phylogeny of Hamatopoecilimon as a step- 
ping stone to further phylogeographic research. We selected two mitochondrial 
protein coding regions — nicotinamide adenine dinucleotide dehydrogenase subu- 
nit 2 (ND2) and cytochrome C oxidase subunit | (COI). These markers, especially 
COI, are widely used in phylogenetic studies of orthopterans (e.g. ALLEGRUCCI et 
al. 2017, KAYA et al. 2015, CHOBANOV et al. 2017). In addition, we test the phylo- 
genetic utility of the mitochondrial control region (CR, A+T- rich region) in Hama- 
topoecilimon. 


Material and Methods 


DNA Sampling and molecular marker amplification 

Material was collected during field trips across the Aegean islands in 2018-2019 
and stored in absolute ethanol or was received from colleages. All species in the 
subgenus Hamatopoecilimon were represented in our sample (Table 1). Total ge- 
nomic DNA was extracted applying salt-ethanol protocol (ALUANABI & MARTINEZ 
1997). A fragment from the 3'end of COI was amplified with forward primer C1- 
J1751- GGAGCTCCAGATATAGCTTTTCC and reverse TL2-N-3014 — TCTAA 
TGCATTAATCTGCCATCTTA (SIMON et al. 1994). ND2 was amplified with for- 
ward TM-J210 — AATTAAGCTAATGGGTTCATACCC, paired with TW-N1284 — 
AYAGCTTTGAARGYTATTAGTTT (SIMON et al. 2006). CR was amplified with SR- 
J14610 — ATAATAGGGTATCTAATCCTAGT and T1-N18 — CTCTATCAARRT 
AAYCCTTT (SIMON et al. 2006). PCR was performed using Thermo Scientific™ 
DreamTaq™ Hot Start PCR Master Mix, Thermo Fisher Scientific Inc., following 
the protocol of the manufacturer. For COI and ND2 we applied thermal cycling 
suggested by CHOBANOV et al. (2017). Cycling protocol from ZHAO et al. (2011) 
was applied to obtain the A+T-rich region, together with a small fragment of 12S 
rDna. PCR product purification and Sanger sequencing was carried out by Macro- 
gen Europe (Macrogen, inc., Amsterdam, The Netherlands) as an external service. 


Sequence preparation and phylogenetic analysis 

Chronograms were visualized and assembled with CodonCode Aligner version 
8.0.2 (CodonCode, Dedham, MA, USA). Sequences were aligned via MUSCLE 
algorithm (EDGAR 2004), implemented in MEGA X (KUMAR et al. 2018). Selection 
of unique haplotypes, scanning protein coding sequences for stop codons and 
substitution saturation tests were carried out in DAMBE ver. 7.2.43 (XIA et al. 2003, 
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XIA 2018). A total of four matrices were created — ND2, COI, CR, and ND2+COl+ 
CR. Best-fit partitioning schemes and models of evolution were selected using 
PartitionFinder ver. 2.1.1 (LANFEAR et al. 2016), run on the CIPRES Science Gate- 
way webserver (MILLER et al. 2010). Results were later used to set phylogenetic 
analysis. Maximum likelihood (ML) analyses were conducted using RAxML ver. 
8.2.12 (STAMATAKIS 2014) (on CIPRES Science Gateway). Node support was ob- 
tained through 1000 bootstrap resampling. Bayesian inference (BI) analyses were 
accomplished in Mr. Bayes 3.2.5 (RONQUIST & HUELSENBECK 2003, RONQUIST et 
al. 2005). Four simulations of Markov chains and 1 X 10° generations were run, 
sampling 1 of every 100 trees. The MCMC parameters were monitored in Tracer 
ver. 1.7.1 (RAMBAUT et al. 2018). The first 25% of trees were discarded as burnin. 


Table 1: Decimal coordinates of localities and GenBank accession numbers of se- 
quences (ND2 - mitochondrial nicotinamide adenine dinucleotide dehydrogen- 
ase gene; COI - mitochondrially encoded cytochrome c oxidase I; CR — mito- 
chondrial control 


[Species [Locality | fe [n [cor Jo 


Poecilimon Chios island, Ns 247022 | 25.955876 | MT416239 | MT416229 | MT416251 
unispinosus | Olimpi 

Poecilimon Chios island, 38.247022 | 25.955876 MT416228 | MT416252 
unispinosus | Olimpi 

Poecilimon Ikaria island 37.63446 26.2887 MT416240 | MT416231 | MT416253 
ikariensis 

Poecilimon Ikaria island 37.63446 26.2887 MT416241 | MT416230 | MT416254 
ikariensis 


Poecilimon Andros island, | 37.8678 24.920 MT416243 | MT416234 | MT416255 
klausger- Apikia 
hardi 


Poecilimon | Andros island, | 37.8678 | 24.920 MT416244 | MT416233 | MT416256 
klausger- Apikia 

hardi 

Poecilimon | Fournoi Island | 37.574 26.4763 MT416242 | MT416232 | MT416257 
hamatus 

Poecilimon | Kos Island, 36.7619 | 26.9601 MN114185 | MN114162 | MT416258 
hamatus Kefalos 

Poecilimon |Naxosisland | 37.0577 ‘| 25.51 MT416245 | MT416235 | MT416259 
paros 

37.0577 | 25.51 MT416246 | MT416236 | MT416260 
paros 


Poecilimon Karpathos 35.5591 27.159 MT416247 | MT416237 | MT416261 
deplanatus island, above 

Volada vill. 
Poecilimon Karpathos 35.5591 27.159 MT416248 MT416262 
deplanatus island, above 

Volada vill. 
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[Species [Locality | e fm fo fo 


Poecilimon Crete island, CAR 37465 23.83878 MT416238 | MT416227 | MT416249 
cretensis Prases 





Poecilimon Crete island, 35.1969 24.2597 MN114198 | MN114173 | MT416250 
cretensis Skaloti 


Results 


Characteristics of the datasets 

The final alignment of ND2 consisted of 957 bp of which 375 variable and 331 
parsimony informative sites. The COI fragment consisted of 1008 bp of which 287 
variable and 258 parsimony informative sites. The final alignment of the control 
region was 976 bp including 15 indels, 295 variable sites and 284 parsimony in- 
formative sites. No stop codons were detected in protein coding sequences. The 
sequence saturation test (XIA et al. 2003) did not show any substantial saturation. 
The concatenated matrix resulted in a dataset of 2941 bp (ND2+COI+CR). 


Phylogenetic analysis 

A total of eight trees were inferred, using different markers and methods. BI and 
ML methods showed agreement on the tree topology (Fig. 1), although ML anal- 
yses show generally lower support (bootstrap). All single gene trees and the con- 
catenated dataset strongly support the early split of P. unispinosus + P. ikariensis 
branch. All trees support the monophyly of the other four species in Hamatopoe- 
cilimon. COI (Fig. 1a) and CR (Fig. 1b), separately, support (P. hamatus+P.depla- 
natus)+(P. paros+P. klausgerhardi) topology, which is in concordance with the 
results from the concatenated matrix (Fig. 1d). All trees support the P. paros+ 
P. klausgerhardi clade. Most gene trees received poor support for a P. hamatus+ 
P. deplanatus clade with the exception of BI analysis of the CR matrix — Bayesian 
posterior probability close to 1 (Fig 1b). The concatenated mitochondrial dataset 
(NAD+COI+CR) (fig. 1d) showed improved support for most nodes. The BI tree, 
inferred from the concatenated dataset, showed full support of all nodes presenting 
the topology (P. unispinosus+P. ikariensis) + ((P.hamatus+P.deplanatus)+(P. pa- 
ros+P. klausgerhardi)). The ND2 tree (Fig. 1c) is the only one showing different 
topology. ND2 suggests P. deplanatus being a sister clade to P. hamatus+(P.pa- 
ros+P. klausgerhardi) with highest support in comparison to the other single ma- 
trices, which corresponds to the nuclear and mitochondrial phylogenies by ULLRICH 
et al. (2010). 
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Figure 1: Mitochondrial phylogenetic trees of subgenus Hamatopoecilimon, based on ML 
and BI analyses of: a) 1008 bp of COI; b) 976 bp of the control region; c) 957 
bp of ND2; d) concatenated dataset (ND2+COI+CR) with 2941 bp. Branch num- 
bers show branch support (ML bootstrap/BI Bayesian posterior probability. 


Discussion 


Phylogenetic performance of molecular markers 

COI gene is a widely used marker for inferring mitochondrial phylogeny, showing 
strong phylogenetic signal on a broad taxonomic scale (SPICER 1995, GOTO & 
KIMURA 2001, REMIGIO & HEBERT 2003, KJER et al. 2014, WANG et al. 2017). It is 
relatively long and easy to amplify due to availability of conserved primer se- 
quences (FOLMER 1994, ZHANG & HEWITT 1997b, SIMON et al. 1994, 2006). ND2 
performed better than other mitochondrial markers in phylogenetic analysis of 
Odonata due to the high proportion of variable and parsimony informative (PI) 
sites, low substitution rate heterogeneity and other performance criteria (CHENG et 
al. 2018). The situation in the genus /sophya, closely related to Poecilimon, is sim- 
ilar as ND2 shows about 20% more variable and PI sites than COI (CHOBANOV et 
al. 2017). Our study confirms that even with small taxon sampling and closely re- 
lated species, ND2 provides more variable and PI sites than COI and overall higher 
branch support. 
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Results from the control region dataset are generally consistent with these from 
protein-coding genes. Both highly variable parts and conserved structural blocks 
are observed in the A+T-rich region, however, in Orthoptera they are not separated 
into distinct domains (compare CLARY & WOLSTENHOLME 1987, ZHANG & HEWITT 
1997a, ZHAO et al. 2011). While variation of the control region above the genus 
level could be extremely high, restraining sequence comparison (MARDULYN et al. 
2003, AMARAL et al. 2016), CR is a useful phylogenetic marker for closely related 
taxa, providing additional information (AMARAL et al. 2016, LI & LIANG 2018), thus, 
with the potential of providing better resolution than COI (VILA & BJORKLUND 2004, 
CHENG et al. 2018). CR was recently used to test speciation hypothesis and detect 
barriers to gene flow between Poecilimon populations (EWELEIT et al. 2015). In our 
study, CR produced a well-resolved phylogeny and provided additional phyloge- 
netic signal in the concatenated mitochondrial dataset. 


Phylogeny and history of lineages 

Mitochondrial phylogeny based on 16S rRNA, tRNA-Val and 12S rRNA fragment 
(ULLRICH et al. 2010) suggested early separation of the northern group (P. unispi- 
nosus+P. ikariensis) which was corroborated by our study. HELLER et al. (2011) 
points out that P. hamatus, P. paros and P. klausgerhardi share the unusual bifur- 
cate shape of the cerci and form the southern group of Hamatopoecilimon. Our 
analyses confirm that P. deplanatus, endemic to the island of Karpathos, shares a 
common ancestor with the southern group. This species, which does not have bifid 
cerci, branches out first on the ND2 tree (Fig. 1c) as well as on the mitochondrial 
tree by ULLRICH et al. (2010), but surprisingly groups with P. hamatus in other 
analyses (Fig. 1a, b, d), although usually with low support. Failure to obtain strong 
node support and unambiguous tree topology is a common pattern, observed in 
phylogenetic studies of terrestrial organisms in the Aegean archipelago. This phe- 
nomenon is explained by simultaneous splits of many ancient lineages, thus, add- 
ing taxa or markers does not improve resolution (compare GANTENBEIN & KEIGHT- 
LEY 2004, PARMAKELIS et al. 2005, KLOSSA-KILIA et al. 2006, ALLEGRUCCI et al. 
2009). 


Phylogeographical implication 

Poecilimon shares many phylogeographic traits with Eupholidoptera — another ge- 
nus that originated on the Aegeid plate (ÇIPLAK 2004, 2010, HELLER et al 2011). 
One of the key geological events that affected the terrestrial lineages in the Aegean 
is the formation of the Mid-Aegean Trench 12-9 Mya (CREUTZBURG 1963, 
DERMITZAKIS & PAPANIKOLAU 1981), which isolated western and eastern lineages 
by sea. The following desiccation of the Mediterranean causing the Messinian sa- 
linity crisis (MSC) 5.9- 5.33 Mya (KRIJGSMAN et al. 2018) reestablished terrestrial 
connections. The only species of Poecilimon on Crete is P. cretensis Werner, 
1903, which is distant from Hamatopoecilimon and diverged later (ULLRICH et al. 
2010). Since Hamatopoecilimon is basal of Poecilimon and there are no older lin- 
eages on Crete or Peloponnese, Poecilimon possibly originated on the east side 
of the Aegean trench. The Aegean islands may thus have been colonized later 
after the appearance of suitable land bridges. Similar scenario is proposed for the 
biogeography of Greek species of Dolichopoda (ALLEGRUCCI et al. 2009). 
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Karpathos was connected to Rhodes and Anatolia during the early Pliocene until 
3-3.5 Mya (DAAMS & VAN DER WEERD 1980). However, some Karpathos lineages 
are directly linked to Anatolia, while distant from these on Rhodes, which may be 
suggesting different age of vicariance events during Miocene (see VEITH & STEIN- 
FARTZ 2004, PARMAKELIS et al. 2006a, b). Situation is similar in Eupholidoptera as 
E. uvarovi (Karabag, 1952) from Karpathos groups with E. akdeniz Unal & 
Naskrecki, 2002 from southeastern Turkey (ÇIPLAK et al. 2009, 2010). This might 
be explained by later re-colonization of Rhodes as during Late Pliocene it was fully 
submerged (POPOV et al. 2004) or at least partly covered with brakish waters 
(DERMITZAKIS 1990). The lineage of P. deplanatus was thus isolated from the 
‘southern’ group during the Middle or Late Pliocene, while the ‘northern’ group 
(P. unispinosus and P. ikariensis) diverged earlier, possibly at the end of the Mio- 
cene. 


Though mitochondrial phylogeny supports simultaneous divergence of P. deplana- 
tus, P. hamatus and P. paros/P. klausgerhardi, morphological traits clearly place 
P. deplanatus in its own lineage. Early divergence of P. deplanatus (before Mid 
Pliocene) combined with polytomy at the tree root of the "southern group" may be 
linked with gene exchange between those lineages prior to the isolation of Kar- 
pathos. Based on genetic distances and phylogeny, dispersal of P. paros+P. klaus- 
gerhardi lineage to the Cyclades may have thus also occurred during the Pliocene 
as suggested for a few lineages (compare POULAKAKIS et al. 2014, KORNILIOS et 
al. 2018). This may explain the lack of Hamatopoecilimon on the Balkan mainland 
and on Crete. Therefore, we have a case of dispersal of the "southern group" 
followed by isolation contrary to the possible vicariance scenario in the case of 
P. ikariensis and P. unispinosus, P. hamatus and P. deplanatus. The split between 
P. paros and P. klausgerhardi may be quite recent as suggested by low genetic 
distances, morphological similarity and lack of acoustic differentiation (HELLER et 
al. 2011). The latter might have happened during Pleistocene as a result of re- 
peated shifts in the sea level and thus connection-disconnection of the Southern 
and Northern Cyclades. 


Less may be speculated about the "northern group". P. ikariensis is present only 
on Ikaria, while P. unispinosus has wider distribution on the Islands of Chios, 
Lesvos and the neighboring mainland Turkey. As Ikaria was recently connected to 
the mainland, its separation may have allowed the "last stand" for P. ikariensis 
escaping from possible competition with P. unispinosus on the mainland (Fig. 2). 


Hamatopoecilimon represents a remnant from an ancient lineage that split com- 
paratively early in the geological history of the Aegean area. This generally agrees 
with the results from the closely related genus /sophya, whose TMRCA was esti- 
mated at 8-9 Mya (CHOBANOV et al. 2017). Data from our study would be useful for 
testing further hypotheses regarding the origin and evolution of Poecilimon and 
other related genera in connection with the paleogeography of the Aegean region. 
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Figure 2: Spatial dimension of the mitochondrial phylogeny of Hamatopoecilimon. 
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